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Abstract 

The green house effect is caused by the gases which absorb infrared ray (IR) emitted by the earth. 
It is worthwhile if we can adjudicate on which gas causes the greenhouse effect in our class. For 
this purpose, one of our authors, Kaneko has designed an educational tool for testing greenhouse 

n 

effect [1[. This system (hereafter abbreviated PAS) is constructed based on photo acoustic effect. 
Without difficulty and high cost, we can build PAS and check the IR absorption of gas. In this paper 
we give the simple theoretical basis for this PAS. The amplitude of sound observed in PAS depends 
on the modulation frequency of IR pulse. Its dependence can be explained by this simple model. 
Further we show the sound amplitude does not depend on the thermal diffusion, which provides 
the accuracy of amplitude as the IR absorption rate of the gas. According to this model, sound 
signal is not the sinusoidal function and it has higher harmonics. The theory and experiment are 
compared in the third harmonics by spectrum analysis. From this apparatus and theory, students 
can study not only the greenhouse effect but also the basics of physics. 
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I. PHOTO ACOUSTIC EFFECT 



The history of photo acoustic (PA) effect is very old. It was found 130 years ago by 
G. Bell and discussed by several famous authors \'^. Recently the Photo Acoustic system 
(PAS) is renewed by Kaneko et al. as an educational tool to show the infrared absorption 

n n 

of greenhouse gases qualitatively [Jj, [3|. PAS is shown in figure 1. The IR source (1) is the 
heat source, 3 Q heating element given 5 V voltage. The IR coming out from source passes 
through the rotating circle with 5 holes (Chopper (2) ), and to be a IR pulse. Then the IR 
pulse enters into the gas cell (3). This IR pulse gives the molecules vibrational energy, and 
its energy is transformed into the heat of gas when no re-radiation occurs. The obtained heat 
pulses give rise to the pressure oscillation and that gives the sound. The sound is transformed 
into electric signal by microphone (5) and its form can be observed on oscilloscope (6). Two 
IR windows attached on gas cell are made of KBr with thickness 2mm, that are not only 
transparent for IR, but also they work to cut the sound from outside. As we can see in figure 
2, the gas cell (3) has cylindrical form with interior radius bo = 10mm, length /q = 100mm, 
and IR passing region has radius b < 10mm. (We change the radius b in experiment) The 
CO2 and H2O (Vapor) provide intense sound signals but such PA signals are not caused 
by the major atmospheric components, N2 and O2. This system clearly demonstrates the 
origin of the greenhouse effect. 

To check the IR absorption of gas directly, it seems plausible to examine the simple mea- 
surement of temperature of gases under the static IR irradiation, however, the temperature 
increase is very small and we cannot obtain the clear result. The reason is the following. 
First, since the pace of temperature increase is too slow, we cannot neglect the effect of 
thermal diffusion to outside. Therefore the gas temperature we observe is not determined 
only by IR absorption of gas. Second, the effect that the gas container is warmed up by 
IR and the gas is warmed up secondary by container might be large. This is related to the 
material of gas container. It is necessary for the gas container that the transmission factor 
for IR to be high enough. Third, the temperature increase is too small and so the usual 
thermometer can not be used. We need to observe the temperature by using a different 
physical variable. 

The PAS solves these three problems and shows the clear IR absorption for greenhouse 
gases. The precise explanation will be given in further section, but we give here some 
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FIG. 1. Schematic diagram of a PA experimental system for IR absorption of gases. 

comments. First, by giving short periodic IR pulse to the gas, we can make the effect of 
thermal diffusion negligible. Second, we use KBr which has very high transparency to IR, as 
gas cell window (4). Third, we observe the temperature variation by the pressure variation, 
that is, the sound. Then we have very sensitive observation without time delay. From these 
reasons we can observe the IR absorption of gas very clearly by PAS. 

This PAS provides opportunities to learn the difference in infrared absorption between 
greenhouse gases and other atmospheric gases for high school level student. However, this 
educational material is still not enough of a challenge to study quantitative analysis. So we 
need some physical quantities which can be observed easily and that can be compared to 
simple theory. For this purpose, we pay an attention to the form and amplitude of sound 
signal of PA. From the experiment, we find two facts for the PA signal. First, though 
IR pulses are given to the cell, the obtained sound signals are similar to sinusoidal waves. 
Second, the amplitude of sound decreases with the increasing modulation frequency of IR 
pulse. In this paper we will explain these two phenomena by using simple phenomenological 
model, which makes quantitative analyses possible in this PAS. Furthermore we show the 
sound amplitude does not depend on the thermal diffusion, which provides the accuracy 
of amplitude to measure the IR absorption rate of the gas. In addition we will present an 
analogy of this PAS to the electric circuit, which helps students for understanding the PA 
phenomenon observed in this educational apparatus. 

II. MODEL 

During the incidental IR beam passes through the gas cell, some of its energy is absorbed 
into the gas molecules as excitation energy of molecular vibrations. Then the molecular 
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FIG. 2. Gas cell has a form of cylinder with inner radius 60 and length /q. The circular IR window 
has radius b with b < bo. 

internal energy is transported into the kinetic energy (center of mass ), and increases the 
gas temperature. We define the effective incidental IR power density S{x,t)[W /m^] in gas 
cell, as is explained in Appendix A. We further define the temperature field 9{x,t) and 
pressure field P{x,t) in the gas cell. Then we obtain two equations for these fields, j4| 

^2n 99 ^ dP 

where we have used the state equation for ideal gas, and the quantities with suffix "0" means 
ambient constant quantities, k, p, and Cp are the heat conductivity of the gas, mass density, 
and the specific heat at constant pressure per unit mass respectively. 

These are simple linear equations, however, it is bit far from intuitive understanding. So 
we hope to obtain the simpler model for this PAS which predicts quantitative experimental 
results. To carry out this idea, we integrate two equations spatially and we adopt the spatial 
mean temperature 6 rather than the local one 6{x,t) as is shown in Appendix B. Then we 
obtain the heat relation. 

where the mean temperature of gas 9, temperature of environment surrounding cavity Oenv 
{6 > 9env), and the heat capacity of gas (isovolumetric heat capacity). M = {2 — 
(6/6o)^}/(8vrZo/«) is the thermal resistance which is explained in Appendix B, where Iq is 
the length of gas cell, Bq is the inner radius of the gas cell, b is the radius of cell window 
(See figure 2). L(t) is the spacial integration of IR power density S, that is, the net energy 
absorption per second as explained in appendix A. 
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FIG. 3. (a) Chopper and front window of the ceh, (b) IR luminosity into the cell 

This equation is nothing but the first law of thermodynamics, which means that the IR 
energy given from outside is distributed into the gas and the generated heat diffuses out of 
cell. Then the gas temperature increases by the difference of these incoming and outgoing 
energy. 

The function L{t) has a form which is similar to the trapezoidal pulse wave as shown 
in figure 3 (b). The reason is given in the following. At t = 0, the edge of front window 
of the cell meets IR coming through a hole region of chopper (figure 4 (a)). After this 
moment, a hole region of chopper partially wrap over the window and L{t) becomes larger 
(figure 4(b)). At t = /i, the window of the cell is just completely covered by a hole region 
of chopper and L{t) becomes maximum (figure 4(c)). In a while fi < t < T/2, L{t) takes 
maximum value (figure 4 (c)-(e)). From t = T/2 to t = T/2 + /i, window of the cell is 
partially covered by blind region of chopper and L(t) goes down (figure 4 (e)-(g)). During 
the time T/2 + fi < t < T, window is covered by blind region of chopper and IR is not 
coming in (figure 4 (g)-(i)). At t = T the situation is the same as Fig. 4 (a). In this way, the 
time period T of IR pulse is corresponding to the angle 0i, and the time /i is corresponding 
to the angle 02 , where 0i is the central angle a set of blind and hole regions of the chopper 
span, and 02 is the central angle the light entrance window of the cell spans (See figure 
3(a)). So we obtain the relation. 



III. SOLUTION 

The initial condition of differential equation ([3]) is given by 
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FIG. 4. The positions of rotating chopper holes and the window of the ceh. This shows the time 
evolution during one time period T. The chopper has 5 holes and so we have 5 IR pulses during 
one rotation of chopper. 

^(0) = ^en.. (5) 

Then the general solution is obtained. 

m = Oen. + L(t')e-(*-*')/-rft', (6) 

lo Jo 

where we used the definition 

r = (7) 

Since the form of L{t) is not simple, we first decompose one to constant part and oscil- 
lation part. 

IR Pulse Constant Part Oscillation Part 




FIG. 5. Decomposition of IR pulse 



L(t) = Lo(l + /(t)), (8) 

where f{t) is the oscillating function with amplitude 1 and period T. By putting this into 
dnj, we have 

m = 0en. + ^Lo{l - e-*/-) + ^e-*/- fity^dt'. (9) 



To perform the last integration 

ffity'/^dt', (10) 

Jo 

we follow the steps below. 

First, f{t) satisfies the relation. 

f{t)=f{t+T), r f{t)dt=Q. (11) 

Jo 

Second, we suppose the condition 

T « r = (12) 

which we will prove later. 

In general we can write t = riT + eT (riiinteger, < e < 1). Then we obtain 

ft , , rnT , , i-nT+eT 

/ f{ty^dt' = / f{ty '^dt' + / f{ty^dt'. (13) 

Jo Jo JnT 



By using e^y^^ = 1 (0 < y < 1) from (dl]) and ([12]), we obtain 

rnT , , 1 _ pnT/T 



/ fity/^dt'= % f{ty'^dt' = n fit')dt' = o, (14) 

Jo 1 — e-^'^ Jo Jo 



and 



/ f{ty^dt' ^ e"^/" / /(t')c^t' = e*/" / f{t')dt'. (15) 

JnT JO JO 

From these two relations we can change eq. as follows. 

m = Oenv + ^i^0(l - e-*/^) + ^ f{t')dt'. (16) 

The physical meaning is the following. The second term in r.h.s. shows the temperature 
change due to constant part. At t << r, temperature grow up linearly with gradient 
Lq^/t, but soon one reaches to the constant value at t » r. This means the energy 
balance between IR incoming energy and outgoing energy by thermal diffusion. 

= Lo. (17) 



On the oscillation part (third part of r.h.s.), because of the rapid oscillation satisfying 
T «T, the heat diffusion to the outside of cavity cannot chase the rapid change. Then the 
integration of luminosity determines the temperature amplitude. This is quite important 
condition to measure the IR absorption rate of gas correctly. 

IV. COMPARED TO EXPERIMENT 

For our approximation to be vahd, we need the condition T « r. First we estimate the 
relaxation time r for our apparatus. The heat capacity ^ is calculated to be 4.06 x 10~^ 
J/K from the sample cell volume of 3.14 x 10~^ m^ and the density of 1.98 kg/m^ and 
isovolumetric specific heat of 28.8 J/(mol • K) for CO2. The thermal resistance M (See 
appendix) is calculated to be 27.4 ~ 54.8 K/W from the cylinder cell's length 1.00 x 10"-^ 
m, the gas's heat conductivity oi k, — 1.45 x 10~^ W/(m • K), and b/bg = ~ l.Then we 
obtain r = 1 ~ 2 s, for b/bo = 1 ~ 0. We work with frequency 1/T = 100 ~ lOOOOHz, that 
is, T = 10^^ ~ 10^^. Then we always have T << r as good approximation. 

A. Form of sound wave 

The IR pulse looks like trapezoidal wave but its integration is very close to a sinusoidal 
function. This is shown in figure 6 theoretically and in figure 7 experimentally. Because the 
rising and falling part of trapezoidal wave is expressed by linear function, and so its inte- 
grated form is parabolic, and the flat part of trapezoidal wave becomes linear function after 
the integration. The combination of these two functions look like the sinusoidal function. Of 
course this is approximately true, and the obtained temperature oscillation includes higher 
harmonics. On this point we will comeback later. 

B. Frequency dependence of sound amplitude 

The normalized oscillating luminosity f{t) can be expanded like 

00 

fit) = ^ a„ sm{nuj{t - /x/2)}, (18) 

n=l 
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FIG. 6. /(t) and its integrated form as temperature oscillation, where amplitude is normalized. 
H/T = 0.24 




FIG. 7. PA signal of CO2 with vertical scale of 100 mV/div and horizontal scale of 10 ms/div. 
where u is the IR modulation frequency. 

Furthermore from the symmetry 

fit + fx/2 + T/4) = fi-t + /i/2 + T/4), (20) 

we have 



a2n = 0. 

By using equation flTB]) . we obtain oscillation part of temperature. 



(21) 



r 00 

mt) = -^T. 7^^. cos{(2n + iMt - ^/2)}. 



u;<^„t^o(2n + l)^^--'^ ' (^^) 
By using the state equation for an ideal gas under the condition of constant volume, we 
obtain pressure oscillation. 



A 



SP{t) = -- E TT?^ cos{a;(2r2 + l)(t - fi/2)}, 



^ ;:ro (2ri + 1 



(23) 



with 



A 



ZRg Lq 



9 



Amplitude [mV] 



200 
180 
160 
140 
120 

100 

80 

60 





















n "XT 


















r 




























































.94 

















































60 90 120 150 180 210 240270 300 

Frequency [Hz] 



FIG. 8. Relation between PA signal amplitude and IR modulation frequency. 

where z is the number of moles of the gas, Rg is the gas constant, and V is the inner volume 
of the cell. The amphtude \5P\ is given as 



A °° 

\SP\ = -5P{t = ^^/2) ^-Y: 



(24) 



The sound amplitude is obviously proportional to the inverse of modulation frequency ou. 

In the experiment, the log-log plot of amplitude versus frequency is given in Figure 8. 
Although the frequency dependence is expressed as f-^-^^ above 100 Hz, which is consistent 
with the theory with 6% error, it is expressed as /~"°'^^ below this frequency. The incon- 
sistency in the low-frequency region can be ascribed to the frequency characteristics of the 
microphone and amplifier in our system. The sensitivity of the microphone and amplifier 
increases with frequency up to about 100 Hz and it is almost constant in the region of 100- 
lOOOOHz. So we can conclude that the amplitude of PA signal is inversely proportional to 
the modulation frequency. 



C. Form of a„ 

The function f{t) is too much complicated and so we can not show the analytic form of 
a„, however, f{t) can be approximated by trapezoidal function. In this case we have simple 
form for a„. 

8 J2n + l)iy. , , 

a2n+\^—n. — nr;2~^^"i n -I' a2n = 0. 71 = 0,1,2 (25) 

'K\2n -|- Ij^i^ z 
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where v is the constant rate factor of IR pulse. 

V = — = /iw. (26) 

Then we have rectangular pulse wave for u = 0, and triangular pulse wave for i/ = tt. In 
practice, the value of u is determined by the relative size of the light entrance window of 
the cell and the chopper. As shown in Figure 3 (a), z/ is expressed as 

u/{2n) = fi/T = 02/01. 

To be exact, the shape of the entrance window gives an secondary effect on the form of 
the IR input power L{t). When the entrance window has a shape of a sector, L{t) has an 
exact trapezoidal form. However, the shape of L{t) becomes a gsmooth trapezoidal wave" 
with round corners for a system having a circular entrance window as our case, which is 
shown in Figure 3. In this case, we can not obtain analytical solutions for the coefficients 
an, so that we need to evaluate them numerically. We will discuss the characteristics of 
sound signals as well as the influence of factor u in both cases in section 5. 

The u dependence of the amplitude can be obtained as follows. Since Lq is proportional 
to the area of cell window vrfo^ as seen in appendix, and 6^ oc 02 ^"^ fixed 0i from figure 
3(a), A has the u dependence as v4(z/) oc z/^. Further we have the relation 



sin{(H!!±l)il} 



to (2n + l)3 ^ 
Then by using ( !24l) and ( !25l) . we obtain the u dependence of the amplitude. 

\6P\ ~ u^{l - u/{2tt))/uj, (27) 

V. HIGHER HARMONICS OF PA SOUND 

The observed PA signal is not an exact sinusoidal wave because it includes higher har- 
monics. In this section we evaluate the contribution of the higher harmonics. Let p„ as 
the n-th harmonics amplitude in PA signal. According to Eq. (123|) . we have for trapezoidal 
approximation, 
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Pi I = ^sin(i./2), (28) 

TTUU 



From this relation we have 



f \-\ I I I 1 + 2 cos z/ I 
^l'^) =1 PsM 1= • (30) 

Let us consider the form of function g[y). giv) is the decreasing function for < z/ < 27r/3 
and increasing function for 27r/3 < < tt. The physical reason is the following. The slope 
of rising edge of trapezoid changes from infinity to smaller value for larger z/, and so the 
higher harmonics are suppressed i^giy) decreases) in a region of small v. But in a region of 
large z/, the peak edge becomes more sharply bent for larger v and the higher harmonics 
enhance {^g[y) increases). This explains the form of dashed line (theoretical value for a sector 
window) in Figure 9. 

The solid line shows the theoretical values for a circular window due to the numerical 
calculation. The qualitative property is the same as the case of sector window, but the 
position of singularity (ps = 0) changes to higher v. 

For the experimental values of spectrum of PA sound with circular window, the first and 
third harmonics appear clearly, but higher ones are buried in the background noise. So 
we can compare only the value of | ps/pi | between theory and experiment. The change 
of V can be controlled by changing the radius of entrance window of cell h with fixed feg 
experimentally. As shown in Figure 9, we have the well coincidence with the experimental 
values and theoretical curve. 5| This shows the validity of our model. For the experiment 
it is necessary to give a comment. Since the IR is not prepared as parallel rays, even the 
shadow part of chopper gains the IR. To avoid such a situation, it is important to put 
chopper very close to the window of the cell, and putting the source of IR as far as possible. 

VI. ANALOGY WITH ELECTRIC CIRCUIT 

To understand a new phenomenon the method of analogy is very convenient. As we have 
already shown, our PAS is very close to the parallel CR circuit. 
From Figure 10, we find the relation. 
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FIG. 9. z/-dependence of harmonics ratio | ps/pi |. The solid and dashed hnes show the theoretical 
values for a circular and a sector window, respectively. The former corresponds to the smoothed 
trapezoidal wave, and the latter to the trapezoidal one. The box, circle, and triangle shaped points 
show the experimental results with circular hole. 



i{t) 




FIG. 10. Parallel CR circuit 



Let us show the correspondence relation between the variables in the electric circuit and 
those in the PAS. 

• total electric current (i) — > IR power (L), 

• electric resistance (R) heat resistance {^), 

• electric capacity (C) — >■ heat capacity (^), 

• electric charge of condenser (g) 

— >■ internal heat energy of gas {U — 't^{9 — 9env)), 

Prom the above relations, we see that the resistive element corresponds to the heat 
resistance against the heat flow into the wall of gas cell, and the condenser corresponds to 
the heat capacity of the gas in the cell. 
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By using the analogy, we obtain the following relation from Eq. (13T|) . which is equivalent 
to Eq. ([3]). Just like the electric circuit where the condenser is charged during a half cycle 
{i{t) > 0) and the condenser is discharged into the resistive element in another half cycle 
{i{t) = 0), the temperature of the gas increases during the gas is irradiated by IR light 
and the temperature decreases by the energy flow during the rest of the irradiation, which 
explains well the IR input and the sound signal shown in Figure 4. 

VII. CONCLUSION 

We have introduced a simple model to deal with the photo acoustic phenomenon of a gas 
enclosed in a cell, whose energy gains and losses can be regarded to have the same relation 
as those in the parallel CR circuit. By using this model we can quantitatively explain the 
following three questions, 

• Why do we observe the sound signal having a sinusoidal like form when trapezoidal 
IR pulses are inputted into the sample gas ? 

• How does the amphtude of the sound signal depend on the modulation frequency of 
the IR input? 

• Why can we neglect the diffusion effect to measure the IR absorption rate in PAS ? 

The analogy to the electric circuit leads us to understanding this phenomenon intuitively. 
Though this model is simple, it is a very convenient way to explain our PAS for students. 
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VIII. APPENDIX 

A. INCIDENTAL IR BEAM 

Incidental IR beam power with unit [W/m^] is expressed like 
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I{x,t)^Ioe-^^{l + m), 

where / is the oscillating function with amplitude 1, and (3 is the absorption rate per unit 
length. Then the power density with unit [W/m^]given to gas is 

-§ = /3V-'^^(l + /(t)). 

Some of their energy is re-emitted and do not contribute to increase the local temperature. 
By using probability of radiation-less transition a, we have an effective incident power 
density, 

We then define the net energy absorption per second [W], 

L{t)^ J S d'x^Lo{l + m), 

with 

Lo = a7r6^/o(l-e"^'«), 
where Iq is the length of gas cell, and nb^ is the area of cell window. 



B. BASIC EQUATION, AND ITS REDUCTION TO MEAN FIELD EQUATION 

Our starting equations are, 

06 dP 
KW^e-poCp— + S = -—, (32) 

and 

Each field is composed of ambient constant and fiuctuation field as, 

e ^eo + 6e{x,t), p = Po + 6P{x,t). 

Then we define the mean field by 
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e{t) = 00 + ^ J^6e{x,t)d^x, (34) 
P{t) =Po + ^ f 6Pix, t)d^x, (35) 



V Jj: 

where S is the inner space of gas cell, and 9S means its boundary hereafter. To obtain the 
equation for these mean fields, we integrate Eqs. ( 132|) and (133|) . 

./^^W.d5-p„Vc,f + L=-Vf. (36) 

and 

On the wall of cell, the normal speed of gas fluid f„ equals zero. 

= Vn = -— fvP-ndt= — VP-n, (38) 

Po J Po^ 

where the linearized Euler equation v = —VP/ p is used, n is the normal unit vector to 
surface of cell, and oscillation with modulation frequency co is supposed. Therefore the first 
term of left hand side of Eq. fl37p equals zero, and so we have after once integration. 

dP_Podl 

dt~ 9, Of ^^^^ 
Then together with Eq. (l36l) . we obtain 

r f)f) 
K Ve-da + L = 't^ (40) 

Jas at 

where the isovolumetric heat capacity = p^Vcp — Nk, with the number of gas particles, 
k the Boltzmann constant. 

Now the outgoing heat flow J is given by 

J= -kVO. 

on the surface of cell inside. We hope to express the heat going out of cell J J ■ da hj 
using 9 but not local function 9(x,t). 

The usual theories for PAS ^ | approximates the coordinate dependence only to the trav- 
eling direction of light x, however, the heat flow goes mainly to the radial direction in our 
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PA sample cell consisting of a metal oblong tube and two KBr windows as shown in Figure 
1. So we approximate our heat flow J going to the radial direction r but not to x direction 
instead. Next we calculate the temperature distribution for static case and consider the 
relation between 9 and temperature gradient at the inner surface of cell. 

Let us consider the static solution of Eq. fl5^ with constant S = Lq/{ti}P'Iq) for r < b 
and S" = for 6 < r < 6o- 

kV^O + S = 0. 

The boundary condition is as follows, drd = at r = 0, and 6 = O^nv at r = where 
6o is the inner radius of cell cylinder. Since the thermal diffusibility of the metal wall of our 
sample cell is high, we treat the temperature at the inside surface of the wall as a constant 
value of Oenv Then the solution of static equation is as follows. First the inner solution 

for r < b. Second the outer solution 

0out{r) = 9{b)-^-^^\^Hr/b), 
ln(6o/o) 

for b < r < bo- To obtain 6{b), additional boundary condition is necessary. 

^ r^ \r=b ^ r\ \r=b ■ 

or or 



This gives solution. 



and 6 is 



eib) = 9,^^, + -^\n{bo/b). 



= —p^ / d{r)rd(j)dr = / rOindr + / rOoutdr) 
irbQ J bo Jo Jb 

.,„„^M2_(vw!}, 

OTTIqK 

Note that from Lq oc 6^, ^ — )■ Oenv as 6 — )■ 0. 

The outgoing heat energy from gas cell per second is the same as net energy absorption 
per second from IR in static case. So we have / J -da = Lq. Therefore we interpret Eq. f l4T|) 

as 

e-eenv = ^ J J-da, (42) 
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with definition of heat resistance 



W^? (43) 



We suppose Eq. (H2!) holds even for the case of non static temperature. Then from Eqs. 
(go]) and g2]), we obtain Eq. (jS]). 
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